Anomalous Roughening of Viscous Fluid Fronts in Spontaneous Imbibition 
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We report experiments on spontaneous imbibition of a viscous fluid by a model porous medium 
in the absence of gravity. The average position of the interface satisfies Washburn's law. Scaling of 
the interface fluctuations suggests a dynamic exponent z ~ 3, indicative of global dynamics driven 
by capillary forces. The complete set of exponents clearly shows that interfaces are not self-affine, 
exhibiting distinct local and global scaling, both for time (/3 = 0.64 ± 0.02, /3* = 0.33 ± 0.03) and 
space (a = 1.94±0.20, ai oc — 0.94±0.10). These values are compatible with an intrinsic anomalous 
scaling scenario. 

PACS numbers: 47.55.Mh, 68.35.Ct, 05.40.-a 



The dynamics of immiscible fluid-fluid displacements 
in porous media has been a subject of much interest in 

last years [IIEHSHIISSI!, both from a fundamen- 
tal point of view, as a dynamical non equilibrium process 
leading to rough interfaces and from a technological 
point of view, in industrial and environmental problems 
such as oil recovery, irrigation, and filtration 0]. The 
process is called spontaneous imbibition when an invad- 
ing fluid that wets preferentially the medium displaces a 
resident fluid at constant external pressure. In contrast, 
forced-flow imbibition takes place when the invading fluid 
gets into the porous medium at constant injection rate. 

The case of spontaneous imbibition has special interest 
from the point of view of dynamic scaling. Spontaneous 
imbibition is dominated by capillarity. In contrast to 
forced-flow imbibition, there is no global mass conser- 
vation l aw, a fact that may have a dramatic effect on 
scaling |lOl Hl| . Spontaneous imbibition is a case with 
slowing-down dynamics in which time scales change con- 
tinuously. Indeed, it is known from a long time |12| that 
the temporal scaling of the mean advancing front obeys 
Washburn's law, (h) ~ t 1 / 2 , but little is known about the 
scaling of fluctuations in this regime. Moreover, most 
experiments in spontaneous imbibition have been per- 
formed with the disordered medium placed vertically , 
a situation where gravity limits the Washburn behavior 
to a short observation time. Experiments with gravity 
have been directed to the study of the pinning process 
and its scaling. Only one experiment with the disordered 
medium placed horizontally has been performed on pa- 
per In this experiment scale invariant roughness is 
found only for highly anisotropic disorder, while loga- 
rithmic roughness, at best, is found in isotropic paper. 
Paper, however, is a medium with uncontrolled disorder, 
and phenomena such as swelling and prewetting might 
change the effect of static capillary forces. 

In this Letter we report on spontaneous imbibition ex- 



periments of a viscous fluid by a model porous medium 
with well controlled disorder, consisting on a horizontal 
Hele-Shaw cell with random variations in gap spacing. 
Our scaling results suggest a dynamic exponent z ~ 3, 
consistent with a nonlocal dynamics driven by capillary 
forces. This exponent has never been observed in experi- 
ments before, despite ha ving been suggested by simple 
models of imbibition 0, Il3| After a detailed scaling 
analysis of the interfacial fluctuations we find different 
global and local exponents, with values that suggest that 
our experiment must be described in the framework of 
intrinsic anomalous scaling [Llj . 

Experimental setup.- The experimental setup 0, Q 
consists on a Hele-Shaw cell of size 190 x 500 mm 2 (L x 
H) made of two glass plates separated a narrow gap of 
thickness b = 0.46 mm. Fluctuations in the gap spacing 
are provided by copper obstacles (d — 0.06 mm high, 
1.50 x 1.50 mm 2 size) that are randomly distributed over 
a fiberglass substrate, without overlap, filling 35% of the 
substrate area (disorder SQ 1.50 in Ref. @). 

In the present set of experiments we use a constant 
pressure device that consists on an oil container of se- 
lectable height H, in the range from —50 to 100 mm, in 
steps of 0.1 mm. A silicone oil (Rhodorsil 47 V) with 
kinematic viscosity v = 50 mm 2 /s, density p — 998 
kg/m 3 , and oil-air surface tension a = 20.7 mN/m at 
room temperature, penetrates into the cell through two 
wide holes drilled in the bottom glass, displacing the air 
initially present. The holes connect to a reservoir in the 
cell, with the disorder plate placed just ahead (Fig-QJ. 

The experimental procedure used to prepare a flat ini- 
tial interface is the following: The reservoir in the cell is 
carefully filled and, using an auxiliary syringe, the oil-air 
interface is pushed gently up to a transverse copper track 
2 mm behind the disorder pattern, and the experiment 
starts. The evolution of the oil-air interface is monitored 
using two CCD cameras. A computer records the images 
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FIG. 1: (a) Sketch of the experimental setup, (b) Schematic 
representation of our Hele-Shaw cell with quenched disorder 
in the bottom plate, (c) Temporal evolution of the oil-air 
interface at H — — 9 mm. 




1000 



2000 3000 
t (S) 



4000 



5000 



FIG. 2: Main plot: Squared average interface position as a 
function of time, for different oil column heights H , showing 
that the imbibition process follows Washburn's law, (h) = 
At 1 / 2 . Inset: Dependence of Washburn's law prefactor A 
(squared) on column height H. 



acquired and stores them for processing. 

Washburn's law.- Different heights of the oil column 
in the container, in the range —9 < H < 15 mm, have 
been explored. One example of the temporal evolution 
of the interface is shown in Fig. ^c) for H = —9 mm. 
As a first result (Fig. [2J we observe that for H > —10 
mm the average interface position obeys Washburn's law, 
(h) = At 1 ' 2 . For H < —10 mm the interfaces recede. By 
Darcy's law (h) ~ Vp ~ H/{h), so that A 2 ~ H, as 
verified in the inset of Fig. [21 

Scaling analysis of the rough fronts.- The statistical 
properties of a one-dimensional interface defined by a 
function h(x, t) are usually described in terms of the fluc- 
tuations of h. More precisely, global fluctuations are mea- 



sured by the global interface width, which for a system of 
total lateral size L scales according to the Family- Vicsek 
(FV) ansatz [3 . However, in experiments one cannot 
usually probe the global width, since system size is fixed 
in standard setups. Instead of global fluctuations, one 
measures the local surface fluctuations, which are more 
easily accessible in the laboratory. Local interface fluc- 
tuations are measured by calculating either the spectral 
density of the interface, S(k,t) = (h(k,t)h(—k,t)), where 
h(k, t) is the Fourier transform of the interface and brack- 
ets indicate averages over independent runs of the experi- 
ment, or the local width, w (l,t) = {([h(x, t) — (/i);] 2 }/) 1 / 2 , 
where (•••); denotes an average over x in windows of size 
I. The local width is also expected to satisfy dynamic 
scaling 



w(l,t) = tPg{l/t 1 /'), 



(1) 



where the most general scaling function [Til Hif is ex- 
pected to behave as 



u a, °" if u < 1 
const, if m > 1 



(2) 



Here (3 and z are the growth and dynamic exponents, 
respectively, and a = (3z is the roughness exponent. So, 
for short scales (I < t 1 ! 7 -) the local width grows as 
with (3* = (a — ai oc )/z, whereas for larger scales the 
growth of the global width, t@ , is recovered. This gen- 
eralization of the scaling form is required to account for 
potentially anomalous roughening (a ^ cv; oc ). Note that 
three independent exponents are necessary now to de- 
fine the universality class. The standard self-affine FV 
scaling [15J is then fully recovered when global and local 
roughness exponents become equal, i.e. a = a; oc , f3* = 0. 
Sometimes a more detailed analysis of the scaling func- 
tions is necessary to completely define the generic class of 
scaling. Then both S{k,t) and w(l,t) must be examined. 
The scaling of the spectrum 



S{k,t) = fc- (2o!+1) s(A;i 1/z ), 
where now the scaling function is 



s(u) 



u 2a+1 if u < 1 



(3) 



(4) 



and a s is the spectral exponent, adds valuable additional 
information [1J| . With the same value of the scaling ex- 
ponents the system can belong to one or another class de- 
pending on the value of the auxiliary exponent a s . For 
instance, with a > 1 the system can be either in the 
superrrough, a s — a, or the intrinsic anomalous class, 
cts = ciioc- This is important because it can inform of 
the existence of different symmetries |l7j . 

Experimental results.- Since slight perturbations of 
the interface at the beginning of the experiment are un- 
avoidable, it is difficult to characterize the scaling of the 



3 




0.01 1 — — — — — 1 

10 100 1000 10000 

t (s) 

FIG. 3: Log-log plot of the interfacial width as a function of 
time, for different window sizes I. The experimental param- 
eters are b — 0.46 mm and H = —9 mm. The straight lines, 
with slopes 0.64 and 0.33, result from a data fit in the scaling 
region. 



interfacial fluctuations at very short times. To mini- 
mize this effect we have always considered the subtracted 
width W(l,t), defined as W(l,t) = [w 2 (l, t) - w 2 {l, 0)] 1 / 2 
i. 

We have focused on measuring the values of the scal- 
ing exponents for an oil column height H = —9 mm. In 
these conditions the imposed (negative) pressure gradi- 
ent nearly compensates the traction due to capillarity, 
so that the interface is driven very slowly. The prefac- 
tor of Washburn's law in this case is A = 1.26 ± 0.15 
mm/s 1 / 2 , so that the average interface velocity (initially 
(h(t = 0)) ~ 0.05 mm/s) decreases by about 84% in the 
experiment. We have performed a total of 6 different ex- 
periments, 2 for each of 3 different disorder realizations. 
The plot of W(l,t) as a function of time and for 9 dif- 
ferent windows of size I is shown in Fig. We obtain 
a power law with different exponents for large and small 
window size, (3 = 0.64 ± 0.02 and (3* = 0.33 ± 0.03. The 
value of this last exponent is very robust and clearly dif- 
ferent from 0. Saturation occurs simultaneously in all 
scales at around 6800 s. Both facts rule out a FV scaling 
and indicate the presence of anomalous scaling. 

The temporal evolution of the power spectrum, pre- 
sented in Fig. ^a), provides the value a s — 0.92 ± 0.11. 
In addition, the observation that the spectra at different 
times are shifted upwards rules out the possibility of su- 
perroughness. The remaining scenario, cti oc = a s , can 
be checked up on the collapse suggested in Eq. The 
result of the best collapse, obtained with (3 — 0.64 (previ- 
ously determined) and z = 3, is shown in Fig.^fb). From 
this collapse we obtain ai oc = 0.94 ± 0.10 as a result of a 
data fit. Having determined the three independent expo- 
nents (3, (3* and a\ oc directly from the experimental data, 
the remaining exponents a — 1.94±0.20 and z = 3.0±0.3 
follow from the scaling relations. We have verified that 
the collapse of the spectrum (not shown in the figure) is 
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FIG. 4: (a) Temporal evolution of the power spectra of the 
interfaces analyzed in Fig. [3] The solid straight lines, with 
slope —2.84, are fits to the data in the scaling regions. The 
vertical dashed line indicates the value q = 0.67 associated 
with the lateral size of the unit copper obstacle, (b) Collapse 
of the interfacial width shown in Fig. |3] using the values j3 = 
0.64 and 2 = 3. The straight line is the result from a data fit. 

also well fitted with the exponents thus obtained. 

To provide more evidences of anomalous scaling, 
we have investigated the possibility of multiscaling, 
C q (l,t) = ([h(x + l,t) - h{x,t)]i) 1 /i ~ l a «. Figure© 
reveals the existence of a lengthscale separating uniform 
scaling from multiscaling. Exponents of the multiscaling 
region are generally not universal, and depend most of 
times on particular choices in the experimental arrange- 
ment. In our case multiscaling is probably due to the 
inhomogeneous local growth in the sharp edge of copper 
obstacles. The exponent of the uniform scaling region is 
universal, and expected to be compatible with cv; oc . This 
is reasonably satisfied by our data, as shown in Fig. [5] 

Discussion.- The deterministic linearized equation 
for the dynamics of a perturbation of the interface, in the 
form of a normal mode of wavenumber k, in spontaneous 
imbibition, reads [H El: h k = -a(b 2 /12^)k 2 \k\h k - 
{h)\k\hk, where /i = pv is the dynamic viscosity of the 
fluid, and (h) is the average interface velocity. Surface 
tension (the term proportional to fc 2 |/c|) is responsible 
for the damping of short wavelength fluctuations, while 
viscous pressure (the term proportional to (/i)|/c|) damps 
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FIG. 5: Multiscaling analysis performed through the compu- 
tation of the q-order height-height correlation function, for 
interfaces at saturation. 

long wavelength fluctuations. 

The two terms, however, cross over at a characteristic 
length £ x = 27r[cr(& 2 /12^)(l/(/i))] 1 / 2 . In spontaneous 
imbibition, (h) ~ i -1 / 2 , and thus the crossover length 
£ x ~ i 1 / 4 . As the interface invades the porous medium, 
Dube et al. 01 have pointed out that the growth of the 
lateral correlation length due to the fc 2 |fc| term, L c ~ 
t 1 / 3 , is effectively dominated by the slower growth of the 
crossover length, £ x ~ t 1 / 4 . One should then observe a 
dynamic exponent z — 4 for length scales larger than £ x , 
a result verified in Ref. . 

In our experiment, however, L c never reaches £ x , and 
thus the observation of a dynamic exponent z ~ 3 is fully 
justified on theoretical grounds. Indeed, (i) at very short 
times L c grows from while, from (h(t — 0)} and the 
values of a, b, and ji, we get £ x ~ 76 mm, which is about 
half the cell width, (ii) At saturation (L c = L) £ x can 
be estimated as follows: L c = (Dt) 1 ^ 3 , where D ~ 10 3 
mm 3 /s is estimated from the inflection in Fig. 0{b), so 
that saturation occurs indeed at a time t s = L 3 / D ~ 
6800 s. The interface velocity at t s can be determined 
now from Washburn's law, and then we get £ x — 190 
mm at saturation, a value that coincides with the size 
of the cell. Thus £ x is only reached at saturation, (iii) 
This result is consistent with the observation that at the 
end of the experiment (h) ~ 0.16 {hit = 0)}, and we get 
again £ x ~ 190 mm. This reflects also that saturation is 
reached shortly before the end of the experiment. 

Conclusions.- Our experiments provide a controlled 
realization of spontaneous imbibition without gravity. 
The interface roughness is a result of the random spatial 
distribution of capillary forces. Traction due to capillar- 
ity is so strong that the external pressure must be nega- 
tive to have a very slow (forward) motion of the rough- 
ening front. Once the control parameter (the pressure 
imposed) is fixed, the front advance verifies Washburn's 
law. 

Our results a ~ 2, a ioc ~ 1, (3 ~ 2/3, j3* ~ 1/3 
and z ~ 3, together with a s ~ ai oc , identify the present 
experiment with the so called intrinsic anomalous scaling 



|14| . This type of scaling is entirely different from the 
usual FV scaling found in the same setup and with the 
same kind of disorder, in forced-flow conditions 6J. This 
substantial difference should be attributed to the absence 
of a global mass conservation law. 

Finally, our results provide the first experimental ev- 
idence of z ~ 3 in spontaneous imbibition. This has 
been made possible by driving the interfaces very slowly, 
thus avoiding the lateral correlation length L c to reach 
the crossover length £ x in the experiment. Preliminary 
experiments in which the interfaces are driven at much 
higher pressure (H — 29 mm) and, consequently, £ x is 
much smaller throughout the whole experiment, cannot 
be properly collapsed with z = 3. 
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